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Chapter 1 

Half-Skyrmion theory for high-temperature superconductivity 

Takao Morinari 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, 

Japan 

We review the half-Skyrmion theory for copper-oxide high-temperature supercon- 
ductivity. In the theory, doped holes create a half-Skyrmion spin texture which 
is characterized by a topological charge. The formation of the half-Skyrmion is 
described in the single hole doped system, and then the half-Skyrmion excita- 
tion spectrum is compared with the angle-resolved photoemission spectroscopy 
results in the undoped system. Multi-half-Skyrmion configurations are studied by 
numerical simulations. We show that half-Skyrmions carry non-vanishing topo- 
logical charge density below a critical hole doping concentration ~ 30% even 
in the absence of antiferromagnetic long-range order. The magnetic structure 
factor exhibits incommensurate peaks in stripe ordered configurations of half- 
Skyrmions and anti-half-Skyrmions. The interaction mediated by half-Skyrmions 
leads to dj.2_j,2-wave superconductivity. We also describe pseudogap behavior 
arising from the excitation spectrum of a composite particle of a half-Skyrmion 
and doped hole. 

1.1. Introduction 

One of the most challenging problems in condensed matter physics is to unveil the 
mechanism of high-temperature superconductivity in the copper oxides. Although 
it has past more than two decades since its discovery,^ no established theory ex- 
ists. The most difficult aspect is to cope with strong electron correlations: The 
undoped system of high-temperature superconductors is an insulator. Contrary to 
conventional band insulators, strong Coulomb repulsion makes the system insulat- 
ing. High-temperature superconductivity occurs by doping holes in such a Mott 
insulator.'^ The pairing symmetry is not conventional s-wave but d2;2_j,2-wave.^ 
It is believed that electron-phonon couplings do not play an essential role in the 
mechanism of high-temperature superconductivity because of the strong Coulomb 
repulsion. Searching for a mechanism based on the strong electron correlation is 
necessary. In this chapter, as a candidate providing such a mechanism the half- 
Skyrmion theory is reviewed. 

The plan of the review is as follows. In Sec. 11.21 we review the structure. 
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electronic properties, and the phase diagram of high-temperature superconductors. 
Then, we describe the half-Skyrmion spin texture in a single hole doped system 
in Sec. 11.31 The half-Skyrmion excitation spectrum is compared with the angle- 
resolved photoemission spectroscopy results in the undoped system. Topological 
character and magnetic properties of multi-half-Skyrmion configurations are de- 
scribed in Sec. 11.41 In Sec. 11.51 we describe a mechanism of dx2-y2-'wave supercon- 
ductivity based on half- Sky rmions. In Sec. 11.61 we describe a pseudogap behavior 
in the half-Skyrmion system. It is shown that the energy dispersion of a composite 
particle of a half-Skyrmion and doped hole leads to an arc-like Fermi surface. 

1.2. Review of high-temperature superconductivity 

Although there are a number of high-temperature superconductors, the essential 
structure is the Cu02 plane. Material differences arise from an insulating layer 
sandwiched by Cu02 planes.^ In the parent compound, nine electrons occupy 3d 
orbitals at each copper ion. In the hole picture, there is one hole at each cop- 
per site. The hole band is half-filled but the system is an insulator because of 
a strong Coulomb repulsion. The system is well described by the spin S = 1/2 
antiferromagnetic Heisenberg model on the square lattice with the superexchange 
interaction J ~ 1500K.'' Experimentally and theoretically it is established that the 
ground state is an antiferromagnetic long-range ordered state."* The structure of 
the Cu02 plane and the arrangement of spins at copper sites in the undoped system 
are schematically shown in Fig. Il.lf a). 

This antiferromagnetic long-range order is rapidly suppressed by hole doping. 
In fact only 2 — 3% doping concentration is enough to kill antiferromagnetic long- 
range order. This critical hole concentration is much lower than the percolation 
limit of ^ 40%. High-temperature superconductivity occurs by introducing about 
0.05 holes per copper ion. A schematic phase diagram is shown in Fig. ILlT b). (In 
this review we focus on the hole doped system and do not discuss the electron doped 
system.) 

In the high-temperature superconductors anomalous behaviors are observed in 
physical quantities for temperatures above the superconducting transition temper- 
ature, Tc.^ The phenomenon is called pseudogap. The Fermi surface observed by 
angle-resolved photoemission spectroscopy (ARPES) in the underdoped regime is a 
truncated, arc- like Fermi surface.^ (See for a review, Ref. 7.) In scanning tunneling 
spectroscopy, a gap like feature appears below the pseudogap temperature T* which 
is higher than Tc-^ For temperatures below T* gap- like behaviors are observed in 
NMR, transport coefficients, and optical conductivity. (See for a review, Ref. 5.) 

In the doped system, because of the strong Coulomb repulsion at each copper 
site^ doped holes occupy oxygen p-orbitals. Hole spins interact with copper site 
spins with strong antiferromagnetic Kondo interaction. Because Kondo interaction 
coupling, Jk, is much larger than J and hole hopping matrix elements, there is 
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Fig. 1.1. (a) Two-dimensional Cu02 plane. Filled circles represent copper ions. Open circles 
represent oxygen ions. Arrows are localized spin moments at each copper site. In the ground 
state, those spins have anti-ferromagnetic long-range order, (b) A schematic phase diagram of the 
high-temperature superconductors. The horizontal axis represents the doped hole concentration 
and the vertical axis represents temperature. AF indicates antiferromagnetic long-range order and 
SC indicates superconductivity. Below a characteristic temperature curve denoted by the dashed 
line, the system shows a pseudogap behavior (denoted by PG). SG and FL indicate spin-glass-like 
state and Fermi liquid state, respectively. 

correlation of forming a singlet pair called the Zhang-Rice singlet.^ 

The strong Jk limit leads to the t-J model. ^ The t-J model has been studied 
extensively. (See for a review, Ref. 10.) In the t-J model, double occupancy is 
projected out. One way to deal with this constraint is to use slave- particle formu- 
lations. Based on the resonating valence bond picture proposed by Anderson, a 
spin-charge separation scenario has been applied to the physics of high-temperature 
superconductors.^" From various physical view points different theories have been 
proposed. There is a view in which incommensurate spin correlations observed 
in neutron scattering are associated with stripe order. (See for a reivew, Ref. 12.) 
Chakravarty et al. proposed d-density wave order as competing order against super- 
conductivity to explain the pseudogap phenomenon. The half-Skyrmion theory 
has some connection with these theories which will be discussed later. 

1.3. Single hole doped system 

The high-temperature superconductors are characterized by a rich phase diagram 
shown in Fig. ll.lT b). Remarkably this phase diagram is essentially controlled by 
a single parameter x, the doped hole concentration. Therefore, to understand the 
physics of high-temperature superconductivity it is necessary to establish how to 
describe doped holes. Here we consider a half-Skyrmion spin texture created by a 
doped hole in an antiferromagnetically correlated spins. 
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As the simplest model we consider the single hole doped system. As stated 
in the previous section, the undoped system is described by the antiferromagnetic 
Heisenberg model on the square lattice, 

H = jJ2Si-Sj, (1.1) 

where the summation is taken over the nearest neighbor sites and the vector Sj 
describes a spin S = 1/2 at site i. Theoretically and experimentally it has been 
established that the ground state is the antiferromagnetic long-range ordered state. ^ 
A convenient description of the state is obtained by introducing Schwingcr bosons"'^^ 
and then describing the long-range ordered state in terms of a Bose-Einstein con- 
densate of those Schwinger bosons. In the Schwinger boson theory the spin Sj is 
represented by 

where the components of the vector cr = {crx,cry,c7z) are Pauli spin matrices. To 
describe the spin S = 1/2, the Schwinger bosons must satisfy the constraint, 
J2 CLCjCT = 1- We introduce a mean field Aij = {CitQi — dlCjt) and introduce a 

Lagrange multiplier Xj to impose the constraint. In the Schwinger boson mean field 

theory, "'^^ we assume uniform values for these quantities as Aij = A and Xj = A. 
The energy dispersion of Schwinger bosons is given by = y/X^ — ^J^A^^ with 
7k = (sinfca; + sinA;j;)/2. In the ground state, A = 2 J A. Bose-Einstein condensation 
occurs^^-i^ at k = (±7r/2, ±7r/2). 

Now we consider a hole introduced in the system. The strong interaction between 
the doped hole spin and copper site spins leads to correlation of forming a Zhang- 
Rice spin singlet® as mentioned in the previous section. If the singlet is formed, 
then the Bosc-Einstcin condensate of the Schwinger bosons is suppressed around 
the doped hole position. Generally if the condensate is suppressed at some point 
in two-dimensional space, then a vortex is formed around that point. The vortex 
solution is found by solving the Gross-Pitaevskii equation. For the Schwinger 
bosons, the vortex turns out to be a half-Skyrmion as shown below. 

For the description of the half-Skyrmion, it is convenient to use the non-linear 
sigma model. ^® Low-energy physics of the antiferromagnetic Heisenberg model is 
well described by the non- linear sigma model, ^® 

« KkBT)-' r r , /o__x 2 

S=^ j dT j dPv 
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(1.2) 



where ps is the spin stiffness and Csw is the antiferromagnetic spin-wave velocity. 
(Hereafter we use units in which h= 1.) The unit vector n represents the staggered 
moment and r is the imaginary time. 

In order to describe the correlation of forming a Zhang- Rice singlet pair between 
doped hole spins and copper site spins, one has to be careful about its description. 
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Obviously forming a static singlet state which is realized in the — > oo limit 
does not work. Because such a simple singlet state contradicts with the rapid 
suppression of antiferromagnetic long-range order by hole doping. If static singlet 
states arc formed, then sites occupied by singlets do not interact with the other 
spins at all. The situation is similar to site dilution, and suppression of magnetic 
long-range order is described by the percolation theory. In other words, considering 
a strongly localized wave function of a doped hole at a copper site is not realistic. 
We need to consider a hole wave function extending over some area so that the 
doped hole spin interacts with the other spins. In fact, numerical diagonalization 
studies of the t-J model show a Skyrmion-like spin texture^° when a hole motion is 
restricted to one plaquette. A similar situation may be realized in Li-doped system 
as discussed in Ref. 21. 

To include the effect of the interaction with the other spins, we formulate the 
correlation of forming a Zhang-Rice singlet in the following way. The spin singlet 
wave function of a copper site spin and a hole spin is described by 



This wave function has the form of superposition of the hole spin- up and copper spin- 
down state, |t)ij \ i)cn spin-down and copper spin-up state, |T)cu- 

In order to include the interaction effect, we consider these states separately and 
construct superposition of them. We assume that the spin state at site j is spin-up 
before the introduction of a doped hole. Under this assumption, the spin-up state 
does not change directions of the neighboring spins. So the system is uniform for 
the staggered spin n. By contrast, the spin-down state at site j creates a Skyrmion 
spin texture characterized by a topological charge. 



where Sxx = £yy = and £xy = —£yx = 1- Following Ref. 22, the Skyrmion solution 
is found by making use of an inequality 



^(IT)hli)cu-|i)hlT)cu)- 



Q = — I (fY£ai3 n (r) • [Sen (r) x S^n (r)] , 





The classical energy satisfies. 




The equality holds if and only if 



daTH ± 



(n X 9/311) = 0. 



(1.3) 



This equation is rewritten in a simple form. If we introduce 
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then Eq. (|1.3p is rewritten as 

{dx T idy ) w = 0. 
This equation is the Caucy-Riemann equation. Noting 



21Xew 2Iniii; Iwl^ — 1 



n 



|u;|' + l' |u;|' + l' \w\' + l 



the solutions satisfying the boundary condition n(rj) — ~z and n(r oo) 
with z the unit vector along the z-axis are the Skyrniion spin texture, 

/ 2rix 2riy r'^ — rj^^ 
n (r) — I — ^ — , — ^ — , — ^ — 

with Q — 1 and the anti-Skyrmion spin texture, 

/ 2rix 2riy — r;^ 

with Q = ~1- 

Now we consider superposition of the uniform state and the Skyrmion state. Al- 
though superposition of the two spin configurations is not the solution of the field 
equation, these solutions suggest that the resulting spin configuration is character- 
ized by a topological charge Q with < |Q| < 1. The value of Q is determined 
by making use of the fact that the antiferromagnetic long-range ordered state is 
described by Bose-Einstein condensation of Schwinger bosons. In the CP^ repre- 
sentation of the non-linear sigma model, ^'^ the U(l) gauge field is introduced by 

= ^*cdi^z^ (1-4) 

where the complex field Zg. is defined through 



In terms of the gauge field a^, the topological charge Q is rewritten as 



n=(4zt)a(;| ). (1.5) 



Q = / ^ {d^ay - dyUx). 



From this expression, we see that a spin configuration with Q corresponds to the flux 
2'kQ in the condensate of the Schwinger bosons. On the other hand, each Schwinger 
boson carries the spin S = 1/2. However, there are no S" = 1/2 excitations in the 
low-energy excitation spectrum. Low-lying excitations are antiferromagnetic spin 
waves which carry the spin one. Therefore, all Schwinger bosons are paired and the 
flux quantum is tt similarly to conventional BCS superconductors.^^ The flux value 
is not arbitrary and Q must be in the form of Q = n/2, with n an integer. Taking 
into account the constraint < |(5| < 1, we may conclude \Q\ = 1/2.^^ The spin 
texture with IQj = 1/2 is called half-Skrymion spin texture because the topological 
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Fig. 1.2. (a) Half-Skyrmion spin texture. Arrows represent the directions of the spin at the copper 
sites. Neighboring spins are almost anti-parallel because of the antiferromagnetic correlations. 
Filled circle at the center denotes the core of the half-Skyrmion. (b) Anti-half-Skyrmion spin 
texture. 

charge is half of the Skyrmion spin texture. The half-Skyrmion spin texture and 
the anti-half-Skyrmion spin texture are schematically shown in Fig. 11.21 

Moving half-Skyrmion spin texture is obtained by applying Lorentz boost on 
the static solution above by making use of the Lorentz invariance of the non-linear 
sigma model. The energy dispersion is 



where Eq ~ ^iips is the half-Skyrmion creation energy. On the square lattice the 
dispersion is transformed into 



Note that the lowest energy states are at (±7r/2, ±7r/2) because the Schwinger 
bosons are gapless at those points in the antiferromagnetic long-range ordered state. 
The half-Skyrmion spin texture are mainly formed by Schwinger bosons around 
those points. 

Now we compare the half-Skyrmion excitation spectrum with the ARPES result 
in the undoped system. The excitation spectrum Eq. (jl.7[) is qualitatively in good 
agreement with excitation spectrum obtained by Wells et al?^ The parameters Cgw 
and Eq are determined from the values for the Heisenberg antiferromagnet. We 
use the renormalized factors Zc — 1.17 and Zp — 0.72, which are estimated from 
quantum Monte Carlo simulations^^'^* and a series expansion technique. Using 
these values, we find that the bandwidth is 1.5 J and Eq = I.IJ. The experimentally 
estimated bandwidth by Wells et al. is 2.2 J. From the fitting of the experimental 
data assuming Eq. (II. 7|) . we find Eq ^ J. 

In the undoped system, anomalously broad line shapes are observed by ARPES. 
Line shape broadening is associated with scattering of excitations by fluctuation 
modes. In the half-Skyrmion theory, half-Skymions couple to spin- wave excitations. 
Describing those spin wave excitations in terms of the gauge field fluctuations line 
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Fig. 1.3. The half-Skyrmion dispersion. Horizontal axis represents positions in the Brillouin 
zone. 

shape broadening is studied by applying a strong coupling analysis. The width of 
the broadening is in good agreement with the experiment. 



1.4. Multi half-Skyrmion configurations 

In the previous section, the single half-Skyrmion has been considered. The most 
important physical quantity carried by the half-Skyrmion is the topological charge. 
The topological charge density, which is defined in the continuum as 

9c (r) ^ (r) • [d^n (r) x dyii (r)] , 
has the following form on the lattice, 

<lc {xj, Vj) = {xj,yj)-[n {xj + 1, yj) - n {xj - 1, yj)] x [n {xj,yj + 1) - n {xj,yj - 1)] 

In the single half-Skyrmion state, qc{xj,yj) has a peak around the half-Skyrmion 
position, and vanishes at infinity. The integration of qd^) leads to the quantized 
value Q = ±1/2. 

If there are many half-Skyrmions, do half-Skyrmions keep topological charge? 
In order to answer this question, we carry out a simple numerical simulation. First, 
we put either an XY-spin-vortex or an anti-XY-spin-vortex randomly. (A similar 
numerical simulation is discussed in Ref. 31.) A multi- XY-spin-vortex configuration 
is defined by 



3 



[x-Xjf + {y-yjf' 



y-V] 



{x - xjf + {y - yj f 



(1.8) 
(1.9) 
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where qj = +1 for XY-spin-vortices and qj — —1 for anti-XY-spin-vortices. A 
doped hole is sitting at each spin- vortex position {xj,yj), and n{xj,yj) — 0. Then, 
a random number which ranges from —0.1 to 0.1 is assigned to the z-component of 
the vector n(xi,yi) except for the nearest neighbor sites {xi ± 1, yi) and {xi,yi ± 1). 
After that, the equihbrium configuration of the vectors n(xi, yi) is obtained by the 
relaxation method. At site {xi,yi), n{xi,yi) is updated by 

n{xe,ye) = ^ [n(xi + l,yi) +n{xe~- l,yi) + n{xi,yi + 1) + n{xi,ye - 1)] . 

The constraint |n(a:£,y£)| = 1 is imposed by taking the normalization after the 
update. This update procedure is carried out over all lattice sites except for the 
hole positions {xj,yj) and its nearest neighbor sites, {xj ± l,yj) and {xj,yj ± 1). 
The resulting converged state is an approximate state for a multi-half-Skyrmion- 
anti-half-Skyrmion configuration. As an example. Fig. 11.41 shows the topological 
charge density distribution at the doping concentration x = 0.107. There are re- 
gions where topological charge density is non-zero. Positive (negative) topologi- 
cal charge density region is associated with half-Skyrmions (anti-half-Skyrmions). 
Non-vanishing distribution patterns are observed for x < Xc with Xc ^ 0.30. For 
X > Xc, half-Skyrmions and anti-half-Skyrmions are heavily overlapped. As a result 
topological charges are cancelled out. Therefore, above Xc the topological nature of 
half-Skyrmions is lost. (For related discussions about the effect of thermally excited 
skyrmions and hole induced skyrmions, see Refs. 32,33.) 




Fig. 1.4. Topological charge density distribution in real space at the doping concentration x = 
0.107 on a 64 X 64 lattice. Positive values are shown in red and negative values are shown in blue. 



Now we discuss magnetic properties of multi-half-Skyrmion configurations. The 
magnetic correlation is investigated by the static magnetic structure factor, 

5(q) = ^f^.,-^)5„,(q). 



a, (3 



Here 
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S{q) is measured by neutron scattering experiments. Introducing, 



S{q) is rewritten as 

5(q) 



1 

N 



|S(q)|^-l|q.S(q)|^ 



If there is antiferromagnetic long-range order, then S{q) has a peak at q = (tt, tt) = 
Q, and the peak height is proportional to the number of lattice sites. From numer- 
ical simulations above, we find that S{q) shows incommensurate peaks at positions 
shifted from q = Q. Furthermore, we find that around x — 0.10 the maximum peak 
height is on the order of the square root of the number of lattice sites. Therefore, 
the magnetic long-range order disappears around that doping concentration. 

The physical origin of the incommensurate peaks is found by studying a regular 
configuration of half-Skyrmions. Taking a vortex-anti- vortex configuration given 
by q[xj,yj) — (—1)'^^^^^, an approximate "antiferromagnetic" configuration of half- 
Skyrmions and anti-half-Skyrmions is obtained by the numerical simulation. Fig- 
ure [TT5Ua) shows the magnetic structure factor of the resulting state at a; = 0.0625 
on a 64 X 64 lattice. The incommensulate peaks are found at [nil ± 2^), 7r(l ± 26)) 
with 5 = 0.125. These peaks are associated with the superlattice formed by half- 
Skyrmions and anti-half-Skyrmions. A stripe order case is shown in ll.Sf b) which is 
obtained by taking a vortex-anti- vortex configuration with q(xj.yj) = (—1)^^ . The 
dominant incommensulate peaks are found at (7r(l ± 26), tt) with 6 = 0.125. 




(a) (b) 

Fig. 1.5. (a) Magnetic structure factor in a roughly antiferromagnetically ordered half-Skyrmions 
and anti-half-Skyrmions. (b) Magnetic structure factor in a stripe-like configuration of half- 
Skyrmions and anti-half-Skyrmions. 



Experimentally neutron scattering experiments show incommensurate peaks at 
q= (7r(l ± 2(5), 7r(l± 2(5)) for x< 0.05 and q= (7r(l ± 25), tt) and q = (tt, 7r(l ± 2(5)) 
for X > 0.05.'^'*"'^^ As shown above, such incommensurate peaks are found in some 
configurations of half-Skyrmions and anti-half-Skyrmions. However, there is a quan- 
titative difference. Experimentally it is found that 6 ~ x.^^'^^ In order to explain 
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this experimental result, it is necessary to consider stripe-like configurations of 
half-Skyrmions and anti-half-Skyrmions. To determine the stable configuration of 
half-Skyrmions and anti-half-Skyrmions, we need to take into account the inter- 
action between half-Skyrmions, which is not included in the numerical simulation 
above. Determination of the stable half-Skyrmion configuration is left for future 
work. 

1.5. Mechanism of d-wave superconductivity 

In the half-Skyrmion theory, doped holes create either a half-Skyrmion or an anti- 
halr-Skyrmion at their positions. A half-Skyrmion or an anti-half-Skyrmion is 
bound to each hole and moves together. There are two ways to formulate the 
effect of half-Skyrmions on the doped holes. One way is to take a doped hole and 
the half-Skyrmion created by the hole as a composite particle. This approach is 
formulated in the next section, and we shall see that the theory leads to pseudogap 
behavior. 

The other way is to include the effect of half-Skyrmions as fields mediating 
interaction between doped holes. By integrating out the half-Skyrmion degrees 
of freedom, we obtain the interaction between doped holes mediated by the half- 
Skyrmions. In this section, we take this approach and show that the interaction 
leads to a dj;2_y2-wave Cooper pairing between the doped holes. ^^''^^ An intuitive 
interpretation is also given about the origin of the pairing interaction based on a 
Berry phase. 

The fact that each doped hole carries a half-Skyrmion is represented by 

V X a = 5-01 (r) V's (r), (1.10) 

where a is the U(l) gauge field in the CP^ model defined by Eq. (|1.4p . The index s 
labels the sign of the topological charge, s = -I- refers to a half-Skyrmion and s = — 
refers to an anti-half-Skyrmion. 

The interaction between the doped hole current and the gauge field is found as 
follows. Doped holes interact with the spins via a strong Kondo coupling, 

3 

Meanwhile, the doped hole motion is described by 

From a perturbative calculation for the tight-binding model describing the Cu02 
plane, we find the Kondo coupling is Jk — leV. Since Jk is larger than t ~ OAeV, 
we first diagonalize the Kondo coupling term. The diagonalization is carried out by 
the following unitary transformation, 

'^3 ^ ^j/i' 
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where 



^3 1 ^j"T 



Under this transformation, the hopping term is 
Extracting the terms inchiding the gauge field a^, we find 

j S=x,y 

(The effect of other terms is discussed in Ref.39.) After Fourier transforming and 
taking the continuum limit, we obtain 

H^^-^^T. E + (1.11) 

k,q 5—x,y 

Here the effective mass m is introduced by f ~ l/2m. 

Now we derive the pairing interaction between the doped holes from Eqs. p.lOp 
and by eliminating the gauge field. In order to fix the gauge, we take the 

Coulomb gauge. In wavevector space, we set a^; (q) = —^a{q) and ay (q) = 
^a{q). From Eq. pTTO)) . we obtain 

Substituting this expression into Eq. (jl.lip . we obtain 

-^i"t — 2^ 72 '5 '^/k's'cr'/k'+q,s'CT' /k+q,s,(T/k,s,(T■ 

Since we are interested in a Cooper pairing, we focus on terms with k + k' + q = 0. 
After symmetrizing the terms with respect to spin and half-Skyrmion indices, we 
obtain 

w ^ kxk' — kyk'^ t I , l\ A A r n 

k^k' s,s',o-,(t' ) 

(1.12) 

This interaction term leads to a pairing of holes as shown below. 
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Following a standard procedure, *° we apply the BCS mean field theory to the 
interaction (|1.12p . The mean field Hamihonian reads, 



H (/k,s,(T /k,-s,-(T /-k,s,(T /-k,-s,-cr 

k,s,cr 

/ c. n a(+) 





^(-)* ^(+) 



V^k,.; 

The mean fields are defined by 





Ck 

-)* 

k, — S. — CT 

(+)*■ 

k, — s.— cr 



k.s,o" k,s.cr 

a'-' a(+' ■ 

K, — s, — a K, — s, — a 



-Ck 









/k,S.cr ^ 
./*k , — s . — cr 
— k,s,cr 
\/-k,-s.-cr / 



A. 



(+) 

k,S,CT 



(1.13) 



(1.14) 



^^■^- = -^"' ^ -j^r^(Ak',-.,-./k'.,.). (1.15) 

k(^k') v«- J 

The mean field (|1.14p is associated with the Cooper pairing between holes with the 
same spin and the same half-Skyrmion index. On the other hand, the mean field 
p.lSp is associated with the Cooper pairing between holes with the opposite spin 
and the opposite half-Skyrmion index. 

If we consider the interaction p.l2p only, then the pairing states described 
by A^'*'j^ and A|^~j^ are degenerate energetically. However, if we include the 
(anti-)half-Skyrmion-(anti-)half-Skyrmion interaction and the half-Skyrmion-anti- 
half-Skyrmion interaction, the Cooper pairing between holes with the opposite spin 
states and the opposite half-Skyrmion indices is favorable. Because the interac- 
tion between (anti-)half-Skyrmions is repulsive and the interaction between half- 
Skyrmions and anti-half-Skyrmions is attractive. So we may focus on the pairing 
correlations described by A|^ ^ ^. 

At zero temperature, the BCS gap equation is 



.(-) 7r«_ ^ k^k'y - kyk'^ 



where i?k is the quasiparticle excitation energy. This gap equation is divided into 
two equations accor 
the gap equation is 



two equations according to the relative sign between s and a. For } s = + ^ 



August 24, 2009 8:35 World Scientific Review Volume - 9.75in x 6.5in 



14 T. Morinari 



For j = 2, the gap equation is 

. (-) _ Ti-i ^ k^k'y - kyk'^ aI,^^ 



Here we consider the case in which both of A|^ ^ and A^ 2 describe the same pairing 
symmetry. Under this condition, we find |A^| = \aI-1\, and so Ei, is given by 

















(1.18) 



A similar gap equation was analyzed in Ref. 41 in the context of the fractional 
quantum Hall systems. Following Ref. 41, we introduce an ansatz, 

A^;! = Afc exp (-z£0k) , (1.19) 

where £ is an integer and 0]^ = tan^^ {ky/k.j;). Substituting this expression into 
Eq. (|1.16p . the integration with respect to the angle di^i is carried out analytically. 
At this procedure, we find that there is no solution for ^ = 0. Therefore, there is 
no s-wave pairing state. Furthermore, the gap equation has the solutions only for 
^ > 0. The gap Afc satisfies the following equation, 

1 f'' „,k'Ak' fk'Y 1 f°° „,k'Ak' f kV 

Afc = — / dk' — + — / dk' — . (1.20) 

2m Jo Ek' \kj 2m Jk Ek- \k' ) ^ ' 

From the asymptotic forms in ^ oo and fc — )■ 0, we assume the following form for 
Aft, 



Ak/sp 



A{k/kF)\{k < kp), 

A{kF/k)\{k>kF), 



where ei;^ is the Fermi energy of holes and kp is the Fermi wave number. The gap 
A is found numerically. The gap A decreases by increasing £. We find A = 0.916 
for £ = 1 and A = 0.406 for £^2. 

The gap equation (|1.17p is analyzed similarly. However, because of the sign 
difference in the interaction the solution has the following form, 

Aj^-^ = Afc exp (z£0k) , (1-21) 

where Afc is the solution of Eq. ()1.20p . As a result, there are two types of Cooper 
pairs with opposite relative angular momentum. 

For p-wave (£ = 1) gap symmetry, the sum of A^^. }_ and A^^. 2 leads to a p^^-wave 
gap which is unstable in the bulk in the absence of symmetry breaking associated 
with spatial anisotropy. Since s-wave gap symmetry is ruled out as mentioned 
above, the lowest energy sate is obtained for dx2^y2--wave gap symmetry. 

The pairing mechanism based on half- Sky rmions is intuitively understood 
(Fig. II. 6p . According to Eq. (ll.lOp . a half-Skrymion, or a gauge flux, is induced 



August 24, 2009 8:35 



World Scientific Review Volume - 9.75in x 6.5in 



morinari't 



Half- Sky rmion theory for high- temperature superconductivity 



15 



around a hole. If another hole passes the gauge flux region at the Fermi veloc- 
ity, a magnetic Lorentz-force-like interaction acts on that hole according to the 
interaction represented by Eq. (jl.lip . A similar pairing interaction is discussed at 
half-filled Landau levels. "'^'''^ In that system the gauge field is the Chern-Simons 
gauge field^^ whose gauge fluxes cancel the external magnetic held fluxes at the 
mean field level. Gauge field fluctuations give rise to a paring interaction between 
flux attached fermions. 



Fig. 1.6. Interaction between doped holes arising from a Berry phase effect associated with the 
gauge flux created by a half-Skyrmion. 



1.6. Pseudogap 

One of the most intriguing phenomena observed in high-temperature superconduc- 
tors is the so-called pseudogap which is observed in various physical quantities.^ 
Here we focus on the pseudogap behavior observed in ARPES. If the hole con- 
centration is lower than the optimum hole concentration at which the transition 
temperature is the maximum, the Fermi surface is not a conventional Fermi surface 
expected from the band theory. Instead, a truncated, an ark-like Fermi surface is 
observed in ARPES. ^ In order to explain this Fermi arc, a standard approach is to 
consider a coupling with some boson modes, such as spin fluctuations or gauge field 
fluctuations associated with phase fluctuations of a mean fleld, expecting self-energy 
effects in the single body quasiparticle Green's function. However, it is not obvious 
that such a conventional analysis leads to qualitatively different physics. 

The half-Skyrmion theory provides a completely different approach. To describe 
the doped hole dynamics which is associated with the spectral function observed by 
ARPES, we need to include the fact that each doped hole carries a half-Skyrmion. 
For that purpose, the direct way is to take a doped hole and the half-Skyrmion 
created by the hole as a composite particle so that the dynamics of the half-Skyrmion 
is included in the doped hole dynamics.**^ 




/ 



August 24, 2009 8:35 World Scientific Review Volume - 9.75in x 6.5in 



16 T. Morinari 

The Hamiltonian describing the hole hopping is 

Ht = -Y.Y. ^y4c,<x + h.c, (1.22) 

where = i for the nearest neighbor sites, — t\ for the next nearest neighbor 
sites, and = t2 for the third nearest neighbor sites. The parameters tijt and 
^2/^ are chosen so that the Fermi surface in the Fermi hquid phase is reproduced.'*^ 
Now let us include the half-Skyrmion dynamics. The half-Skyrmion dispersion 
is described by the Hamiltonian, 



where the k-summation is taken over the reduced Brillouin zone, ± < tt, 
Cckcr = (cko- + Ck+Q,<T)/\/2 and Cokcr = (ckcr - Ck+Q,o-)/\/2, and 

Kk = —V [(cos fca; + cos ky) + i (cos fca; — cos /Cj^)] . 

The half-Skyrmion dispersion is given by ±|Kk|. This dispersion corresponds to 
Eq. (|1.7p with £'0 = and Csw — v. The half-Skyrmion creation energy Eq is zero 
because it vanishes in the absence of the antiferromagnetic long-range order. The 
spin-wave velocity for the doped system is denoted by v which is different from Cs„ in 
the undoped system. Here we use the same creation operators and the annihilation 
operators for doped holes and half-Skyrmions. Because a doped hole and the half- 
Skyrmion carried by the hole is taken as a composite particle. Note that it is not 
necessary to distinguish between half-Skyrmions and anti-half-Skyrmions. Because 
their excitation spectra are the same. 

The dispersion energy of the composite particle is calculated from H = Ht+ i/hs 

as, 

-(-) 



Kk 



with e'^^ — (ek±ek+Q)/2 and ek — —2t (cos kx + cos ky) — 4ti cos k^ cos ky — 
2t2 {cos2kx + cos2ky). The spectral function is calculated following a standard 
procedure.*^ The imaginary time Green's function for up-spin is defined as 

GkT (r) = -(r.CkT {r)cl^ (0)) , 

where Ckt(T) = e'^'^^^^^-'ckte^'^'-^"^^-' with /i the chemical potential and N the 
number operator and T^- is the imaginary time (r) ordering operator. 
The Matsubara Green's function is obtained by 

GkT (i^n) = r^^^ dre'^-^G^^ (r) , 







where aj„ = 7r(2n + l)fcsT (n — 0, ±1, ±2, ...) is the fermion Matsubara frequensy. 
By the analytic continuation, iujn uj + iS, with S an infinitesimal positive number. 
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the retarded Green's function is obtained. Thus, the spectral function is 



1 



All (cj) = ImGki {l^ + iS) 



(1 + Ck) S (u 



(+) 



(1 - Ck) S {io 



^l) , (1.23) 



where, 



Ck 



(2t + v) (cos + cos ky) 



\J {2t + w)^ (cos kx + cos kyf' + v"^ (cos kx — cos ky 



In numerical computations, the parameter 5 above is taken to be finite so that 
the i5-function in the right hand side of Eq. (|1.23p is replaced by the Lorentz function. 
Figure fTTTT a) shows Ak(w = 0) with 5/t — 0.10. The resulting Fermi surface is arc- 
like because the factors 1 ± Ck suppress the intensity in part of the Brillouin zone. 
This arc-like Fermi surface is consistent with the ARPES results. 

However, there is some deviation from the experiment around the ends of the 
arc.'*^ Basically the Fermi arc follows the underlying Fermi surface which appears 
at high-temperature above the characteristic pseudogap temperature, T* . The de- 
viation is suppressed by including the effect of the short-range antiferromagnetic 
correlation. ^^'^^ The effect is included by taking average over the wave vector 
change in the (incommensurate) antiferromagnetic correlation. The calculation is 
similar to that described in Rcf. 49. Figure fTTTy b) shows how the spectral intensity 
is modified at the antiferromagnetic correlation length £^af = 10. 




k,J IT 




(b) 

Fig. 1.7. (a) The intensity plot of the spectral function A\^{w = 0) in a quadrant of the Brillouin 
zone. The doping concentration is a; = 0.10. The other parameters are ti/t = —0.25, t2/t = 0.10, 
and v/t = 1.0. (b) The intensity plot of the spectral function A-i^(u: = 0) with including the effect 
of the short-range antiferromagnetic correlation. The antiferromagnetic correlation length is taken 
as E^AF = 10. 



Now we comment on the similarity to the d-density wave theory. -'^^ The calcu- 
lation for the spectral function A\^{uj) is almost identical to that in the d-density 
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wave theory.^^ However, the physical interpretation of v is different. In d-density 
wave theory, v is associated with d-density wave long-range order. But here v is 
associated with the half-Skyrmions dynamics. Although the origin is different, the 
half-Skyrmion theory and the d-dcnsity wave theory may share some results about 
the pseudogap phenomenon because the Hamiltonian is almost the same. 



1.7. Summary 

In this chapter, we have reviewed the half-Skyrmion theory for high-temperature 
superconductivity. We have mainly focused on four topics. First, we have discussed 
the single hole doped system. The correlation of forming a singlet state between 
a doped hole spin and a copper spin has been investigated including the fact that 
the hole wave function extends over a space, and affected by the other copper site 
spins. A half-Skyrmion spin texture created by a doped hole has been described. 
It is shown that the half-Skyrmion excitation spectrum is in good agreement with 
the ARPES results in the undoped system. 

The most important aspect of the half-Skyrmion is that the doped hole carries a 
topological charge which is represented by a gauge flux in the CP^ representation. 
This property is in stark contrast to conventional spin polaron pictures. Because in 
that case interaction clouds arising from spin correlations are not characterized by 
a topological charge or a gauge flux. 

Second, we have considered multi-half-Skymion conflgurations. In multi-half- 
Skyrmion configurations, antifcrromagnctic long-range order is suppressed around 
the doping concentration, x — 0.10. The topological property of half-Skymions is 
kept for X < 0.30 as shown by numerical simulations. After suppression of antifer- 
romagnetic long-range order, the spin correlation becomes incommensurate. Nu- 
merical simulations suggest that the origin of the incommensurate spin correlation 
is associated with stripe conflgurations of half-Skyrmions and anti-half-Skyrmions. 

Third, we have discussed a mechanism of d^2_j^2-wave superconductivity. The 
gauge flux created by half-Skyrmions induces the interaction between doped holes. 
The interaction leads to a d^2_y2-wave superconducting state of doped holes. The 
origin of the attractive interaction is a Lorentz force acting on a hole moving in a 
gauge flux created by another hole. 

Finally, we have discussed a pseudogap phenomenon. We have considered a 
composite particle of a hole and half-Skyrmion. The pseudogap is associated with 
the excitation spectrum of the composite particle. 

Although several aspects of the half-Skyrmion theory for high-temperature su- 
perconductivity are described in this review, we need further studies to establish 
the theory. In particular it is necessary to show the half-Skyrmion formation in 
the single hole doped system in a more convincing way to provide a sound starting 
point. 
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